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ABSTRACT 

We present long slit spectrophotometry considering the presence of thermal 
inhomogeneities (t 2 ) of two HII regions in the Small Magellanic Cloud (SMC): 
NGC 456 and NGC 460. Physical conditions and chemical abundances were 
determined for three positions in NGC 456 and one position in NGC 460, first 
under the assumption of uniform temperature and then allowing for the possibil- 
ity of thermal inhomogeneities. We determined t 2 values based on three different 
methods: i) by comparing the temperature derived using oxygen forbidden lines 
with the temperature derived using helium recombination lines, ii) by compar- 
ing the abundances derived from oxygen forbidden lines with those derived from 
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oxygen recombination lines, and iii) by comparing the abundances derived from 
ultraviolet carbon forbidden lines with those derived from optical carbon recom- 
bination lines. The first two methods averaged t 2 =0.067 ± 0.013 for NGC 456 
and t 2 =0.036 ± 0.027 for NGC 460. These values of t 2 imply that when gaseous 
abundances are determined with collisionally excited lines they are underesti- 
mated by a factor of nearly 2. From these objects and others in the literature, 
we find that in order to account for thermal inhomogeneities and dust depletion, 
the O/H ratio in low metallicity HII regions should be corrected by 0.25-0.45 
dex depending on the thermal structure of the nebula, or by 0.35 dex if such 
information is not available. 

Subject headings: H II regions: abundances — galaxies: H II regions H II regions: 
thermal inhomogeneities — H II regions: individual (NGC 456, NGC 460) 



1. Introduction 



The Magellanic Clouds play a fundamental role in astrophysics because they are the 
closest laboratory we have to probe the extragalactic distance scale and the theories about 
chemical evolution of stars and galaxies. They are close enough to be able to identify indi- 
vidual objects and isolate them from their environment. In the case of chemical abundance 
studies of H II regions, it is possible to perform observations that avoid bright stars, Wolf- 
Rayet stars, and supernova remnants; this can translate in extremely good equivalent widths 
in emission, EW em (H/3) > 250 A. 

There have be en many studies of the chemical composition of HII regions in the Mag- 
ellanic Clouds (e.g . 
Caplan et al.lll996l ; 



Peimbert fc Torres- Peimbertlll976t iPagel et al.lll978t IWesterlund 111990 



Peimbert et al.l 120001 ) . It has been found that the chemical composition 



within each cloud is uniform, and that the SMC has lower metallicity (approximately 0.40 
dex) than the Large Magellanic Cloud (LMC); making the SMC the best object to study 
low metallicity environments. In particular, H II regions of the SMC provide the means to 
study low metallicity star- forming neighborhoods and N GC 456 is one of the brightest H II 
regions of the SMC (IPeimbert &: Torres-Peimbertlll976f ). In this study the number of lines 
presented for NGC 456 and NGC 460 is considerably larger than in previous works. 

Almost all studies of HII regions use the direct method to determine gaseous abun- 
dances, and implicitly assume that this value represents the chemical abundances of the ISM. 
The direct method uses collisionally excited lines (CELs) to obtain chemical abundances, 
and presuppose a homogeneous temperature structure. When abundances are calculated us- 
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ing recombination lines (RLs), they are approximately factors of 2 higher than abundances 
obtained with CELs; this is now known as the abundance discrepancy factor (ADF) prob- 
lem. Moreover, the depletion into dust grains can account for a large fraction of the oxygen 
atoms. 



Esteban et al.l (119981 ) measured the depletion into dust in the Orion nebula and sug- 
gested that a 0.08 dex correction to the O/H ratio should be included when considering 
HII regions. More recent studies of the Orion nebula sugg est that the correction is higher 
(IMesa-Delgado et al.l l2009bl ; ISimon-Dfaz fe Stasihskal 1201 ll ) while recent studies of deple- 
tion of Mg, Si, and Fe suggest that the correction depends on the O/H ratio of the specific 
H II region (IPeimbert fe Peimbertll2010l ). reaching a fraction of up to 28% for high metallicity 
regions. 



Peimbertl ( 119671 ) proposed to take into account t 2 when using CELs to determine abun- 



dances in modeling H II regions, and found that spatial thermal fluctuations of approximately 
20% are capable of producing errors of factors of 2 to 3 on the determination of chemical 
abundances. The additional complications to calculate abundances using t 2 , along with the 
fact that there is no tailored physical model to fit temperature inhomogeneities, have made 
its reception within the community poor. 

Since the idea of Peimbert's temperature inhomogeneities was introduced, several ex- 
planations about where such temperature inhomogeneities could be coming from have been 
proposed: shock waves, shadowed regions, advancing ionization fronts, multiple ionizing 
sources, X- rays, magnetic reconnection , and inhomogeneous chem ical composition, among 
others (e.g. IPeimbert fe Peimbertl 120061 ; iTsamis fe Pequignotll2005l . and references therein). 
However, most of these explanations are capable of either producing only a fraction the 
t 2 value found in observations or, under fine-tuned circumstances, the total observed t 2 val- 



ues. 



For most objects, two explanations are capable of producing the observed t 2 values: 
a) the observed t 2 values come from a variety of mechanisms in a chemically homogeneous 
medium, where each mechanism is only responsible for a fraction of the total; and b) the 
presence of chemical inhomogeneities within the photoionized region. As of now, no complete 
self-consistent model has been proposed for either explanation: models for explanation (a) 
have not been able to explain where all the required energy to keep these thermal inhomo- 
geneities comes from, while models for explanation (b) have not been able to describe the 
exact characteristics and evolution of the proposed high metallicity inclusions. As for abun- 
dance determinations: option (a) would indicate t hat RLs give the correct a bundances, or 
that one could use CELs and the formalism of t 2 of IPeimbert fe Costerol (119691 ) to determine 
abundances; and models of option (b) imply that the chemical abundance is intermediate 
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between those obt ained from CELs adopting t 2 =0.00, and RLs adopting a chemically inhomo- 



geneous medium (ITsamis fc Pequignot 



20051 ). For either option, the abundances determined 



using CELs with t 2 =0.00 need a systematic correction. 

Additional evidence supporting chemically homogeneous models with large t 2 values 
implied by the ADF problem comes from: the theoretical vs. observational AY/ AO ratio, 
abundances of F and G young stars of the solar vicinity vs. abundances of HII regions, 
oxygen abundances of O and B stars vs. abundances of HII regions, protosolar abun- 
dan ces vs. HII region abunda nces, considering the effect of Galactic chemical evolution 



see 



Peimbert &: Peimbertll201ll and references therein). 



In Sections [2] and [3] the observations and reduction procedure are described. In Section 
H] temperatures and densities are derived from four and three different line intensity ratios, 
respectively; also in this section we derive the mean square temperature inhomogeneities, t 2 , 
in three different ways: i) from the comparison of the abundances of O II lines to those of 
[O III] lines, ii) using He I lines, and iii) from the comparison of the abundance of C II lines 
to those of [C III] A 1907 + C III] A 1909 lines. In Section Owe determine ionic abundances 
with three different methods: i) based on RLs, ii) based on CELs, and iii) correcting CELs 
for t 2 \ while the total abundances derived from these methods are presented in Section |6j 
In Section [7] we present the direct method and its shortcomings, the typical correction to 
the O/H ratio for low metallicity HII regions, and the implications for the determinations 
of extragalactic chemical abundances. Finally, the discussion and conclusions are presented 
in Sections E] and [9J respectively. 



2. Observations 

The observations were obtained during the night of September 10 th of 2002 with the 
Focal Reducer Low Dispersion Spectrograph 1, FORS1, at the Very Large Telescope, VLT, 
Melipal in Cerro Paranal, Chile. We used three grism settings: GRIS-600B+12, GRIS- 
600R+14 with filter GG435, and GRIS-300V with filter GG375 (see Table 1). 

The slit was oriented almost east-west (position angle 82.5°) to observe the brightest 
regions of NGC 456 (a=l h 13 m 50.8 s , <5=-73°18'03.2") which we called NGC 456 Position 2. 
The linear atmospheric dispersion corrector, LACD, was used to keep the same observed 
region within the slit regardless of the air mass value. The slit length was 410" and the slit 
width was set to 0.51", this setting was chosen to have the resolution to deblend the [O II] 
line A 3726 from A 3729 when using the GRIS-600B+12. Three observations were made with 
each grism configuration to be able to identify and remove cosmic rays. 
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The aperture extractions were made for a length of 5.6" for NGC 456 Position 2; NGC 
456 Position 1 is centered 11.4" east of Position 2 and has a length of 17.2", and NGC 456 
Position 3 is centered 19.4" west of Position 2 and has a length of 22.4" (note there is a 3.2" 
gap to avoid a bright star between positions 2 and 3, see Figured]). The eastern part of 
the slit passed through the outer parts of NGC 460; we defined NGC 460 Position 1 as the 
region centered 199" east of NGC 456 Position 2, with a length of 40" . 

The spectra were reduced using IRAF0 reduction packages, following the standard pro- 
cedure of bias subtraction, aperture extraction, flat fielding, wavelength calibration and flux 
calibration. For fl ux calibration the standa rd stars LTT 2415, LTT 7389, LTT 7987 and 
EG 21 were used (Ifiamuy et al.lll992l . Il994l ). The observed spectra are presented in Figure 
121 a zoom-in of the region around A4600 of the spectrum of NGC 456-2 is presented in Figure 
E] to show the quality of the data. 



3. Line Intensities and reddening correction 

The procedure of line measurement was done using the splot task of the IRAF package. 
Line intensities were measured by integrating the flux in the line between two limits over 
the local continuum estimated by eye. There were a few cases where line blending was 
found. In this situation we used a multiple Gaussian profile procedure. The error in the 
flux calibration has been estimated to be 1%; the contribution to the errors due to the noise 
was estimated from the continuum. The final adopted errors were estimated using standard 
error procedures. 

The reddening correction, C(H/3), and the underlying absorption, EWnh sCH.0), were fit 



ted simultaneously to the theoretical ratios. We adopted the extincti on law of ISeaton ldl979h. 



and t he underlying absorption ratios obtained from the work by I Gonzalez- Delgado et al. 



( 1199a ) 



(1999) and from the theoretical ratios determined with the program INTRAT by lStorey fc Hummer 



The Balmer and helium emission lines were corrected for underlying absorption, and 
the equivalent widths in absorption adopted are presented in Table |2j The extinction law 
values used are presented in Tables |3] and HI The physical conditions used for INTRAT were 
T e =10,000 K and n e =100 cm -3 , this is only a first guess, but since the hydrogen lines are 
nearly independent from temperature and density no further corrections were necessary. 

The emission line intensities of all positions in NGC 456 and NGC 460 are presented in 



1 IRAF is distributed by NOAO, which is operated by AURA, under cooperative agreement with NSF. 
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Tables [3] and HI respectively. Tabled column (1) presents the adopted laboratory wavelength, 
A, column (2) the identification for each line, column (3) presents the extinction law value 
used for each line ( ISeaton Ill979l ). /(A). Columns (4-6) present the data from NGC 456-1, 
showing: the observed flux relative to H/3, F(X); the flux corrected for reddening relative to 
H/3, /(A); and the percentage error associated with those intensities, respectively. Columns 
(7-9) and (10-12) present the same data as columns (4-6) but for positions 2 and 3, 
respectively. In Table H] the columns are the same as the first 6 columns of Table |3] but for 
Position 1 of NGC 460. The EW abs (R(3), C(H/3), EW(E/3), as well as the F(E/3) and I(/3) 
for each position in both objects are also presented in these tables. 



4. Physical Conditions 
4.1. Temperatures and densities 



The sources for the atomic data of CELs in the used IRAF are presented in Table |5] - 
for studies of nebulae with hi gher densities , we re c ommend to use mo r e rece n t Q + collisional 
streng ths s uch as t hose from I Wang et al.l ( 12004 ) ; iMontenegro et al.l (120061 ) ; iPradhan et al. 
( 120061 ). and iTayal I ( 120071 ). The temperatures and densities in Table [6] were determined using 
the line intensities presented in Tables [3] and HI These determinations were done with the 
temden task in IRAF, which models populations in five-, six-, or eight-level ions to derive 
the physical conditions. 



4.2. Temperature Inhomogeneities 

Instead of assuming a homogeneous temperature throughout the objects we took into 
account inho mogeneitie s in the temperature structure, which are described by the formalism 



developed by iPeimbertl (119671 ) . To derive the ionic abundance ratios we used the average 
temperature, T , and the mean square temperature inhomogeneities, t 2 , defined as follows: 

J N e (r)N ion {r)dV 

{ ] ~ T 2 1 N e (r)N ion (r)dV ' W 
where N e and N ion are the electron and the ion densities, respectively, of the line of sight, 
and V is the observed volume. 
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Two independent temperatures are required to infer To and t 2 : one that weights pref- 
erentially the high-t emperature reg ions and another that weights preferentially the low- 
temperature regions (jPeimbertl 1 1 96 71 ) . The temperatu r es tha t weight preferentially the high 
temperature region were defined as in lPeimbert et al.l (120021 ): 



- 4363/5007 



1 + 



t 2 /91300 



(3) 



2 V T 

for the ++ region; for the + region we used a similar equation for the [O II] AA 3727 and 
7325 lines with an energy corresponding to 97800 K. 



Following the literature (jPeimbertl 120031 ; iPeimbert et al.l 120051 . and references therein) , 
we can determine a t 2 value from the helium lines. The helioIO program determines a 
a temperature that weights preferentially the low temperature regions from the He I lines. 
This temperature can be determined with the helioIO program, which uses an analysis 
involving a maximum likelihood method (see Section |5TT|) . which finds the best simultaneous 
fit to T(He I), n e , t 3889 , and He + /H + . This T(He I) can be combined with the temper- 
atures derived from the oxygen CELs to derive t 2 and T . In positions 1 and 3 of NGC 
456, we simply adopted the t 2 value obtained from the He I lines: t 2 =0.035±0.032 and 
t 2 =0.040±0.040, respectively. For NGC 456-2 and NGC 460-1 the values we obtained with 
the helioIO program are t 2 =0.053±0.017 and t 2 =0.032±0.032, respectively. 

We used the temperature derived from the ratio of the multiplet 1 of O II lin es to the 
CELs of [O III] as given by equations [8] — [12] in the work of (IPeimbert et al.ll2005l ). 



T, 



(OII/[OIII]) 



T 



4651/5007 



/i(T ,t 2 



(4) 



With the oxygen RLs in NGC 456-2 we obtained that T =10,120 K and t 2 =0.083±0.019 
using the equations above. For the case of NGC 460-1, T =ll,370 K and t 2 =0.041±0.027 also 
with oxygen RLs. In NGC 456 - 2 the oxygen is about 80% twice ionized (see Section |5~3|) . 
hence the value of t 2 (He + ) ~ t 2 (0 ++ ). NGC 460 has approximately 60% of its oxygen twice 
ionized, but since we only have one other t 2 determination to compare it with, we will also 
assume that t 2 (He + ) « t 2 (0++). 

The combined values of thermal inhomogeneities from the helioIO program and from 
the oxygen RLs resulted in t 2 =0.067±0.013 for NGC 456 Position 2 and t 2 =0.036±0.027 for 
NGC 460 Position 1. 

We were able to obtain the C ++ ionic abundance from both UV CELs and from optical 
RLs (see Section [5]) . We found that the abundances are not consistent with a homogeneous 
temperature structure. A possible value for t 2 (C ++ ) in NGC 456 is 0.08±0.04; this value 
represents the t 2 of the UV observed region. 
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The t 2 determinations in this paper, with the exception of NGC 456-2, lie within typical 
values of ther mal inhomogeneities measured for Galactic H II regions which range between 
0.03 and 0.04 iGarcia-Roias fc Estebanll2007h . The t 2 determination of NGC 456 Position 2 
lies within typical values for extragalactic HII regions, which range between 0.03 and 0.11. 
This is probably because in regions that are farther away it is impossible to isolate "simple" 
volumes and one observes the sum of many physical processes: shadows, high density knots, 
shock waves, ionization fronts, photoionization, etc. 



5. Ionic Abundances 



5.1. Helium 

We used the helioIO program t o obtain the ionic abu ndance of He + , presented in 
Table [7J This program is described in iPeimbert et al.l (120111 ) and it uses a maximum like- 



lihood method to search for the set of parameters (T , t 2 , n e , r 38 89, T(He I), and He + /H + ) 
that produce the optimal simultaneous fit to all the observed helium lines as well as the 
measured T e [0 III] and T e [0 II]. The effective recombinat ion coefficients for the H I and 



He I used were those g i ven b y IStorey fc Hummer! (119951 ) for H and by Benjamin et al. 
( 19991) and iPorter et al.l ([20071 ) for He. The col lisional contribution was estimated from 



Sawey fc Berringtonl (119931 ) and lKingdon fc Ferlandl (|1995l). The o p tical depth effects in the 
triplets were estimated with calculations made by lBenjamin et al.l (120021 ) . 



Correction for underlying stellar absorption is also important for helium lines, in a 
similar way as is for hydrogen lines. Hence, to cor rect for underlying He I absor ption for 
lines with A <5000 A we used values deter mined by Gonzalez-Delgado et al.l ( 1999l ). and for 



redder lines we used the same values as in Peimbert et al 



(120051 ). 



Between 9 and 12 He I lines were used as input in the helioIO program to determine 
the t 2 value for each observed position in both NGC 456 and NGC 460. 

The t 2 derived from He I lines can be used to redetermine abundances in the high- 
temperature regions considering temperature inhomogeneities, due to the similar ionization 
potentials of He + and ++ . Since the He II line A 4686 was not detected, the fraction of 
He ++ /He+ is smaller than lxKT 3 for NGC 456 Position 1, 4xl0" 4 for NGC 456 Position 
2, and 3xl0~ 3 for NGC 456 Position 3 and NGC 460 Position 1. 
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5.2. C and O from Recombination Lines 



The ionic abundance for C ++ presented in Table was obtained from the measurement 
of the A 4267 line of C II. The effective recombination coefficient used was that determined 



by iDavey et al.l (120001 ) for Case A and with T =10,000 K. Although we did measure the A 
4267 line in positions 1 and 3 of NGC 456 and in Position 1 of NGC 460, the errors were 
too large to provide any information other than an upper limit. 

The effective recombination coefficients for determin ing the ionic abundan ce of oxygen 
from t he lines of the multiplet 1 of O II were taken from iPeimbert et al.l (119931 ) and IStorey 
(11994J ) assuming Case B for T e =10,000 K and N e =100 cm -3 . The multiplet consists of 
eight lines each of which de pends on the electron density even though the sum of their 



inten sities, /(sum), does not (IRuiz et al.ll2003l ; iBastin fc Storey! 120061 ; IPeimbert fc Peimbert 



20101 ). Since these lines are very faint, it is often necessary to estimate the unobserved and/or 
blended lines. 

The measurements of the O II lines in NGC 456-2 and NGC 460-1 are presented in 
Tables [3] and HI the abundances derived from these lines are also presented in Table [3 The 
errors in the measurements of the oxygen RLs in positions 1 and 3 of NGC 456 were too 
large to provide useful information. For the first two mentioned cases, we were able to detect 
four of the eight lines in the multiplet. Due to the spectral resolution of the observations, 
those lines were blended into two pairs: AA4639+42 and AA4649+51. 

Measurements of RLs of heavy elements are not common in long-slit spectrophotometry. 
Contamination from blends of other lines of other ions could b e a concern. To check this, we 



looke d into echelle observations of HII regions with good S/N (|Peimbertll2003l ; lEsteban et al. 



2004 ; Garcia-Rojas et al. 2004 ; Garcia- Roj as fc Esteban 20071 ); the only other lines that 
appear in this region of the spectra are: N III AA 4641 and N II A 4643. The contribution 
of the sum of these two N lines compared to the sum of the four O lines observed is about 
10% for objects with metallicity and N/O close to solar, and less than 4% for 30 Doradus 
that has N/O about three times less than solar and slightly higher than the SMC; given our 
S/N, we will ignore it (see Figure [3]). 



5.3. Ionic Abundances from Collisionally Excited Lines 



The ionic carbon abundance from CEL s was determin e d wit h UV spectra of NGC 456 
obtained from the on-line data presented by iBonatto et al.l (119951 ) . We measured the [C III] 
A1 907 C I I I1A190 9 lines using the O III] AA1660+66 lines with the UV C(H/3) presented 
(119791 ). and the atomic data contained in the task ionic of IRAF v2.15.1. The 



in 



Scat on 
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value of the C ++ /H + abundance is presented in Table [8] with and without considering thermal 
inhomogeneit ies . 

With the exception of carbon, all other ionic abundances presented in Table [S] were 
determined using the task abund of IRAF considering only the low and medium-ionization 
zones, which correspond to the low and high-ionization zones of the present work. For posi- 
tions 1, 2, and 3 of NGC 456 we used T low =13,500 K, 12,400 K, and 12,600 K, respectively; 
and T high = 12,650 K, 12,165 K, and 12,300 K, respectively 

Ionic abundances considering thermal inhomogeneities, t 2 ^0.000, were obt ained using 

the tra ditional determinations, 1 2 =0.00 0, corrected by the formalism presented by lPeimbert fc Costero 
dl969h . see also IPeimbert et al.l f|2004h . 



6. Total Abundances 

In general we observed that all positions in NGC 456 presented a high O ionization 
degree implying that there is no substantial amount of He . The amount of He + is similar in 
all positions of both NGC 456 and NGC 460, which may indicate that the ionization degree 
is also similar in Position 1 of NGC 460. 

In the case of C and N, total gaseous abundances were determined with the following 
equations: 

^-ICFfC) ^ ^ (5) 
N(R)- lLHL} A(H+) (5) 

and 



where we used the ionization correction factor (ICF) for C with respect to O of iGarnett et al. 



(1995), and multiplied it by the fraction of C ++ /0 ++ in NGC 456-2 following IPeimbert 



(120031 ) ; this value amounted to 1.27. The N total gaseou s abundance was determined adopt- 



ing the predicted ICF(N) of IPeimbert & Costero! fll969h . N{O)/N{0 



The total gaseous abundances for CI, S, and Ar were determined with the following 
equations: 

N(S) T ^ , A(S+) + A(S++) 

iV(H) =ICF(S) N(jF) ' () 
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and 



iV(Ar) 



ICF(Ar 



N(Ai++) + N(At +++ ) 
iV(H+) ' 



(9) 



wher e the corresponding IC Fs were estimated from iGarnettl (119891) for C I, in th e case of S 
from iPeimbert et al.l (120051 ). and for Ar the ICF was taken from lLiu et al.l ( 120001 ) . For NGC 
456-2 these values amounted to 1.26, 1.36, and 1.30, respectively; and 1.26, 1.12, and 1.71 
for NGC 460-1. 



Following IPeimbert fc Costerol (119691 ) . total gaseous abundances of O and Ne were cal- 
culated with the equations: 



N{0) N(0+) + N(0 + +) 



N(R) 



and 



N(Ne) _ JV(0+) + N(0 ++ ) N(Ne ++ ) 

X 



iV(H) 



iV(0" 



iV(H+) 



(10) 



11) 



Total gaseous abundances for all available elements are presented in Tables 191 and [TD1 for 
NGC 456 positions 1, 2, and 3, and in Table ITTl for NG C 460 Position 1. All these ab undances 
are corrected by depletion of O in dust according to IPeimbert &: Peimbertl ( 120101 ) . Within 
the errors, abundances in all three positions of NGC 456 are very similar. Observations with 
higher S/N of NGC 456-1, -3 and of NGC 460-1 would allow to measure a more accurate 
value of t 2 . 



According to IPeimbert fc Peimbertl (120101 ). a correction of 0.10 dex should be made to 
the total O abundances due to depletion of oxygen in dust. The adopted total abundance 
values for all available elements are presented in Tables [TJ] and H3J For comparison, we also 
present protosolar abundances and abundances o f the Orion nebula, N GC 6822, and NGC 
346, which is the brightest HII region of the SMC (IPeimbert et al.ll2000l ). In addition to this, 
we obtained an abundance determination for the ISM of the solar vicinity from the Galactic 
gradient derived from H II regions, that amo unted to 12+log(0/H) = 8.81 once corrected by 
0.12 for dust depletion (jEsteban et al.l 120051 1 . 



7. A Systematic Correction to the Direct Method 

7.1. The Direct Method and its Shortcomings 

The ADF is defined as the ratio of abundances determined using RLs vs. abun- 
dances obtained with CELs. ADFs are now well established for both Planetary Nebulae 
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and HII regions. RLs have a temperature dependence given approximately by I oc T e 1 , 
thus are brighter at lower temperatures; whereas CELs have a much stronger dependence, 
I oc exp(kT e / AE)T e , and thus brighter at higher temperatures. 



Several studies ( 


Peimbert et al. 


1993 




1995 




Tsamis et al. 


2003 




2004 




Peimbert et al. 


2005 




Garcia-Roias & Esteban 


2007; 


Peimbert et al. 


2007 




Esteban et al. 


2009, and refer- 



ences therein) have shown that ADF typical values are from 1.5 to 3 for most HII regions, 
and from 1.5 to 5 for most Planetary Nebulae. This is clear evidence that the traditional 
or direct method of abundance determinations - involving the assumption of homogeneous 
temperature determined from the e mission line ratio 4363 / (4959+5007) A of [O III] - needs 
to be corrected. Since the work of iPeimbertl (119671 ) and iPeimbert fc Costerd (119691 ) . tem- 
perature inhomogeneities have been p resented as a viable solution to the so called ADF 
problem. IPeimbert &: Peimbertl ( 1201 ll ) discuss several methods that have been used to ob- 
tain t 2 values based on: i) the comparison between the temperatures derived from the ratio 
of the Balmer and Paschen continua to the Balmer line intensities and the temperatures 
derived from CELs; ii) the comparison of the C ++ abundances derived using [C Illj+C III] 
with abundances obtained from C II; hi) the comparison of ++ abundances derived with 
[O III] with abundances obtained using O II; iv) the comparison of the temperature that can 
be derived from He I lines to the temperatures derived from CELs; v) a high s patial resolu- 



tion m ap of 1.5 xlO 6 columnar temperatures in the Orion nebula determined by lO'Dell et al. 



( 120031 ); and vi) the comparison of the different ways to calibrate Pagel's method. In all cases, 
the t 2 values are consistent with each other and the ADF problem goes away. 

In the case of photoionization models, the intensity of nebular lines (such as A 3727 
[O II], or A 6584 [N II]) can be adjusted in most models of both Planetary Nebulae and 
HII regions; however, auroral lines (such as A 4363 [O III]) are in general not adjusted in 
those photoionization models. 



7.2. Typical Correction for Low Metallicity HII regions 

We find that O/H abundances in HII regions determined with the direct method are 
underestimated and have to be corrected by factors of 0.25-0.45 dex. The correction we 
propose is due to two distinct and important physical processes: thermal inhomogeneities 
and depletion of oxygen into dust grains. 

When allowing thermal inhomogeneities to exist in the abundance determination pro- 
cess, we find that the derived oxygen abundances in HII regions increase by 0.15-0.35 dex, 
depending on the particular characteristics of the thermal structure of each photoionized 
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region. 



According to iPeimbert fe Peimbertl (120101 ) , abundances of photoionized regions require 



a correction of 0.09 - 0.10 dex for low metallicity H II regions due to depletion of oxygen into 
dust. Therefore the total correction amounts to about 0.35 dex in the total O abundance. 

Many objects have t 2 determined with uncertainties larger than about 0.030, and fre- 
quently these determinations are consistent with t 2 =0.000; the fact that many such objects 
exist does not mean that the t 2 =0.000 or that there are many objects with negligible tem- 
perature inhomogeneities. Since all the objects where the determinations have uncertainties 
smaller than 0.015, are not consistent with a homogeneous temperature, it means that most 
objects have meaningful temperature inhomogeneities. To ignore the presence of these in- 
homogeneities because the data quality is not good enough, produces a very large bias in 
the study of HII regions in general. For this reasons we recommend to use this average 
correction of 0.35 dex in the total O abundance, and to explicitly present it so that, if the 
recommended correction improves with time, it is clear what to do. 



7.3. Implications for the Determinations of Extragalactic Chemical 

Abundances 



Studies of abundances in H II regio ns can also be used to c heck o ur u nderstanding of the 
chemi cal evolution of our own Galaxy. IPeimbert fe Peimbertl (120111 ) and ICarigi fe Peimbert 
( 120111 ) find that careful determination of chemical abundances in HII regions (once cor- 
rected for depletion of oxygen into dust and for temperature inhomogeneities), are consistent 
with: i) the protosolar abundances after correc ting them for the Galactic chemical evolution 



( Asplund et al 



2009 



Simon-Diaz fe Stasihska 



Carigi fe Peimbertll201lf ). ii) young F and G sta rs of the solar vicinity 



(Bc nsbv fe Feltz ing 2006 ), and iii) O and B stars of the Orion region (jPrzybilla et al.l 12008 



201 lh . 



Since the afore mentioned methods of determining abundances for stars and HII re- 
gions are entirely independent, the consistency in the results allows to individually validate 
each method. In particular it implies that a model with temperature inhomogeneities is 
a better approximation to reality than the direct method, that considers a homogeneous 
temperature structure. 

Almost all studies of H II regions in the literature are made assuming homogeneous, or 
nearly homogeneous, temperature throughout the whole object, hence, the abundances de- 
rived from these works are underestimated by a factor of approximately 2. The study of abun- 
dances in low metallicity H II regions is important to set limits in the models of various areas 
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of astrophysics such as the radial abundance gra dients in spiral galaxies ( IVila-Costas fc Edmunds 



1992 



Zaritskv et al. 



2011 



An et al 



2008 



994 



Esteban et all 120051). the chemical evolution of starburst galaxies 



Peimbert Peimbertl 



lation (e.g. 



Tremonti et al. 



2004; 



Perez-Montero et all 120081 : Bordalo fc Telled l201ll ; ICarigi fc Peimbert 



2011 ), the mass-luminosity relation and the mass-metallicity re- 



Kewlev & Ellison! 120081 : Mannucci et al.ll2010l : iThuan et al. 



2010), the initial mass func t ion, star formation rate and y ield determinations (e.g. lCarigi et al. 



20061 : IPilvurin et al. 



2007; 



dance (IPeimbert et al. 



Peimbert fc Peimbertl 120111). and the primordial helium abun- 



2007; 



Izotov 



2010; 



Peimbert et al.l |2010| . and references therein), 



among others. Consequently, it is of great importance for these studies not to use abun- 
dances derived from the direct method, and to use abundances obtained with a method that 
either allows the presence of thermal inhomogeneities or is not affected by them (e.g. CELs 
corrected by temperature inhomogeneities or RLs). 



8. Discussion: the Thermal Structure of Gaseous Nebulae 



In photoionization models, the temperature is nearly homogeneous because the processes 
considered for both heating and cooling are nearly proportional to the density squared. 
Strong thermal inhomogeneities require processes that generate either heating or cooling 
that are not proportional to the density squared. 

Obs ervationally, the A DF problem emerges when a homogeneous temperature is as- 
sumed. lO'Dell et al.l (120031 ) presented a high spatial resolution map of the columnar elec- 
tron temperatures of the Trapez ium with which they fo und a t 2 value from different lines of 
sight within the nebula, whereas iTsamis fc Walsh! (120111 ). studying components and different 
velocities, found that each line of sight has thermal inhomogeneities. 

All available explanations in the literature for the presence of the ADF require the pres- 
ence of thermal inhomogeneities. There are two families of mechanisms that produce tern- 



i.e. chemical inhomogeneities ( r 


Ibrres-Peimbert et al. 


1990l: 


Tenorio-Taelel 


1996; 


Liu et al. 


2000; 


Tsamis & Pequignot 


2005 


Ercolano et al. 


2007; 


Stasihska et al.l 


2007 


), and ii) an in- 



homogeneous heating, e .g. Shockwaves, advancing p hotoionization fronts, shadowed regions, 
magnetic reconnection (IPeimbert fc Peimbertl 120061 . and references therein). There is a small 
fraction of Planetary Nebulae and H II regions that are chemically inhomogeneous; these ob- 
jects will not be discussed in this paper. Models that are chemically inhomogeneous predict 
abundances that are intermediate between those derived with the direct method and from 
RLs. For most objects an inhomogeneous heating mechanism is in better agreement with 
abundances derived from stellar objects (see section 173]) . and with the thermal structure of 
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the Orion nebula (see section I7T|) . Moreover, t 2 values derived with independent methods 
are consistent with each other within the errors. Specifically, the fact that the t 2 values from 
He and O are consistent, implies that the material in the nebulae is well mixed since these 
elements have different formation histories ( ICarigi fc Peimbertll201ll ). 



Typical photoionization model values of t 2 are about one order of magnitude smaller than 
the temperature inhomogeneities determined from observations. With near photoionized 
objects such as t he Orion nebula, ex tremely detailed studies can be made to measure precise 
values of t 2 : e.g. lO'Dell et al.l (120031 ) presented a high spatial resolution map of the columnar 
electron temperatures of the Trapezium with which they found th at shadowed regions can 
sign ificantly contr i bute to the observed thermal inhomogeneities; iTsamis fc Walsh! (1201 if ) 
and Tsamis et al. (|2011 ) were able to isolate individual high densi ty re gions and found that 
CELs u nderestimate abundances; and iMesa-Delgado et al.l (j2009al ) and iMesa-Delgado et al. 
( 12009bl ) studied the head of HH 202 and found that the ionization front has a t 2 value 
significantly higher than that of the nebula. 

Photoionization models include an energy injection proportional to the density squared; 
this should be a good approximation if the energy were dominated by photoionization. How- 
ever, the measurements of the filling factor (which is defined as the occupied volume over 
the total volume) argue against such an energy budget model. Typical filling factor values 
for HII regions are in the 0.01-0.10 range. These values are showing that HII regions are 
pretty much empty, but they are not showing the origin of the energy that is keeping the 
material apart. This is strong evidence that the traditional energy budget model to study 
HII regions is incomplete. The fact that it is not known what is feeding the filling factor 
in HII regions does not mean that abundance results considering t 2 >0.00 are incorrect, it 
means that models have to be improved to better understand the physics and astrophysics of 
these objects. Temperature inhomogeneities are showing that photoionization models have 
to be corrected to account for the additional physics that has not been included. 

Chemical abundances determined using CELs strongly depend on the thermal structure 
of the nebula, whereas abundances obtained with RLs do not. Even though we know many 
processes capable of generating temperature inhomogeneities, we do not know, quantitatively, 
the relative importance of each one of them. We can estimate the magnitude and the 
effect that these thermal inhomogeneities have on the chemical abundances. Although it is 
interesting to study the origin of such thermal inhomogeneities, it is beyond the scope of 
this paper to do so. In order to determine abundances we need to know the value of t 2 and 
not what is producing it. 
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9. Conclusions 

A detailed analysis was performed for the SMC HII regions NGC 456 and NGC 460. 
This analysis involves abundance determinations considering thermal inhomogeneities (t 2 ^0.00); 
for the purpose of comparison with other authors, we also determined abundances with the 
direct or traditional method, w hich assumes a homog eneous temperature structure (t 2 =0.00). 



The abundances of NGC 346 flPeimbert et al.ll2000h . NGC 456, and NGC 460, are similar 



once corrected for dust depletion. The O/H abundances of the SMC HII regions are 0.58 
dex lower than those of the ISM of the solar vicinity, and 0.46 dex lower than the solar 
photosphere. This indicates that the Small Magellanic Cloud is chemically less evolved than 
the Milky Way and so is the best laboratory to study in detail the ISM of low metallicity 
galaxies. 

Determinations of oxygen abundances in this work were done with three methods: i) 
using RLs, ii) using the direct method for CELs, and iii) using CELs relaxing the homo- 
geneous temperature assumption. The first two methods are not compatible. Temperature 
inhomogeneities of some sort are required to reconcile observations of RLs and CELs. The 
ADF only becomes a problem if one demands a homogeneous thermal structure from the 
start. Our results confirm that the ADF problem can be solved by taking into account the 
presence of thermal inhomogeneities in a chemically homogeneous medium. 

Abundances of H II regions are used to constrain models of many areas of astrophysics, 
therefore, we need to derive abundances from gaseous nebulae with the highest precision 
possible. A crucial aspect of the available observations is that they require a complex thermal 
structure; consequently, thermal inhomogeneities must be taken into account to determine 
abundances. 



There are many ways to obtain a t 2 value (IPeimbert fc Peimbertll201ll ). the most com- 
monly used is based on the comparison of abundances derived with [O III] lines with abun- 
dances derived using O II lines. In this paper we also obtained t 2 values from He I lines 
and from the comparison of abundances derived using IUE UV C III] and [C III] lines with 
abundances derived from optical C II lines. For both NGC 456 and NGC 460, we combined 
t 2 values of ++ and He + to obtain an average t 2 . The t 2 value for C ++ in NGC 456 was not 
included in such average of t 2 because optical and UV observations do not come from the 
same place in the H II region; nonetheless, this t 2 (C ++ ) is relevant because it is not consistent 
with a homogeneous temperature structure in the H II region. 

To obtain accurate values of t 2 , good S/N is required for the faint lines of He I, C II, 
and O II. He I and O II for the same volume are consistent within the errors. Since each 
t 2 value comes from a different element and they have different formation histories, this 
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is evidence that the material within the nebula is well mixed; otherwise those t 2 values of 
different elements would differ from one another. 

The homogeneous temperature assumption of the direct method should only be taken 
as a first approximation because, in addition to photoionization, other energy sources have 
to be taken into account. Temperature inhomogeneities measure the overall importance of 
the many mechanisms that inject energy into the system. We find that O/H abundances 
of low metallicity H II regions are systematically underestimated and have to be corrected 
by a factor of 0.25-0.45 dex. We suggest that, if no additional information is available, an 
increase of 0.35 dex on the O/H ratio should be used; this correction includes temperature 
inhomogeneities and oxygen depletion into dust. 

We are grateful to the anonymous referee for the careful reading of the manuscript and 
several excellent suggestions. MAPG and AP received partial support from UNAM (grant 
PAPIIT 112911), and AP and MP received partial support from CONACyT (grant 129753). 
MTR received partial support from Proyecto Basal PB06 (CATA) and Centro de Astrofsica 
FONDAP. 
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Table 1. Journal of Observations 



Grism Filter A (A) Resolution (A/AA) Exp. Time (s) 



GRIS-600B+12 3450-5900 1300 3x720 

GRIS-600R+14 GG435 5250-7450 1700 3x600 

GRIS-300V GG375 3850-8800 700 3x120 
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Table 2. Absolute Equivalent Widths 



Hydrogen Helium 

ljme EW abs (H/3) ■ blIle EW abs (H/3) 



Ha a 


0.90 


3820 a 


0.108 


H 7 a 


1.05 


4026 a 


0.100 


H<5 a 


1.08 


4388 a 


0.084 


H7 a 


0.99 


4471 a 


0.179 


H ga,6 


0.93 


4922 a 


0.107 


H9 a 


0.78 


4009 d 


0.100 


H10 c 


0.67 


5016 d 


0.114 


Hll c 


0.54 


5048 d 


0.115 


H12 c 


0.42 


5876 d 


0.138 


H13 c 


0.35 


6678 d 


0.082 


H14 c 


0.30 


7281 d 


0.032 


H15 c 


0.25 






H16 c 


0.20 






H17 c 


0.16 






H18 c 


0.13 






H19 c 


0.10 






H20 c 


0.08 






H21 c 


0.06 






H22 c 


0.05 







Gonzalez-Delgado et al.l ( 119991 ). 



b Note that He 1(3889) is included with 
H8. 

c Extrapolations based o n VLT echelle o b- 
servations of 30 Doradus (jPeimbertl 120031 ) . 

d M. Cervino, private conversation. 
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Table 3. NCG 456: Line Intensities for Positions 1, 2, and 3 



Position 1 Position 2 Position 3 



A (A) 


Id. 


/(A) 


F(A) 


7(A) 


% err. 


F(X) 


7(A) 


% err. 


F(X) 


7(A) 


% err 


3587 


He I 


0.272 


0.56 


0.64 


18 


0.59 


0.66 


10 


0.93 


1.02 


28 


3614 


He I 


0.271 


0.95 


1.08 


14 


1.09 


1.22 


7 


0.91 


0.99 


29 


3634 


He I 


0.270 


0.89 


1.01 


15 


1.23 


1.37 


7 








3676 


H22 


0.269 








0.35 


0.40 


13 








3679 


H21 


0.268 








0.43 


0.48 


11 








3683 


H20 


0.267 








0.38 


0.43 


12 


0.68 


0.76 


33 


3687 


H19 


0.266 








0.54 


0.61 


10 


0.64 


0.72 


34 


3692 


H18 


0.265 








0.63 


0.71 


9 


0.55 


0.62 


37 


3695 


H17 


0.264 








0.76 


0.99 


9 


0.99 


1.35 


28 


3704 


H16 


0.262 








1.24 


1.56 


7 


1.18 


1.62 


25 


3712 


H15 


0.260 








1.22 


1.58 


7 


1.55 


2.11 


22 


3722 


H14+[S III] 


0.257 


1.54 


1.99 


11 


1.60 


2.05 


6 


1.46 


2.11 


23 


3726 


[O II] 


0.256 


65.02 


73.16 


2 


58.79 


65.23 


1 


61.13 


66.59 


4 


3729 


[O II] 


0.255 


86.89 


97.72 


2 


73.92 


81.98 


1 


83.15 


90.54 


3 


3734 


H13 


0.254 


2.13 


2.69 


9 


1.70 


2.20 


6 


1.69 


2.44 


21 


3750 


H12 


0.250 


2.38 


3.00 


9 


2.48 


3.11 


5 


2.09 


2.97 


19 


3770 


Hll 


0.245 


3.14 


3.97 


8 


3.24 


4.06 


4 


3.05 


4.20 


16 


3798 


H10 


0.238 


4.11 


5.14 


7 


4.31 


5.34 


4 


3.95 


5.39 


14 


3820 


He I 


0.233 


0.75 


0.92 


16 


0.80 


0.97 


8 


1.65 


1.91 


21 


3836 


H9 


0.229 


5.97 


7.29 


6 


5.98 


7.26 


3 


5.48 


7.19 


12 


3867 


[Ne III] 


0.222 


17.91 


19.82 


3 


22.23 


24.30 


2 


22.82 


24.55 


6 


3889 


H 8+Hc I 


0.218 


16.71 


19.25 


3 


16.57 


18.93 


2 


17.26 


20.04 


7 


3967 


[Ne III] 


0.201 


20.25 


22.68 


3 


22.15 


24.01 


2 


21.25 


22.68 


6 


4009 


He I 


0.193 








0.17 


0.69 


18 








4026 


He I 


0.190 


1.37 


1.58 


12 


1.60 


1.90 


6 


1.71 


1.97 


21 


4067 


[S II] 


0.182 


1.01 


1.10 


14 


0.93 


0.99 


8 


0.71 


0.76 


33 


4076 


r n TTi 

[S 11] 


0.181 


0.32 


0.35 


24 


0.30 


0.32 


14 








4102 


H<5 


0.176 


23.81 


26.65 


3 


23.82 


26.45 


2 


23.62 


26.53 


6 


4125 


[Fe II] 


0.172 








0.30 


0.32 


14 








4144 


He I 


0.169 








0.30 


0.42 


14 








4169 


Oil 


0.164 








0.18 


0.19 


18 








4185 


Oil 


0.161 








0.12 


0.13 


21 








4267 


ClI 


0.143 


<0.30 


<0.30 




0.04 


0.04 


40 


<0.30 


<0.30 




4340 


H 7 


0.128 


45.03 


48.25 


2 


44.52 


47.51 


1 


44.29 


47.28 


4 


4363 


[QUI] 


0.122 


4.12 


4.34 


7 


4.36 


4.57 


4 


4.22 


4.38 


13 


4388 


He I 


0.116 


0.58 


0.66 


18 


0.42 


0.50 


12 








4471 


He I 


0.094 


3.77 


4.13 


7 


3.72 


3.97 


4 


3.89 


4.20 


14 


4563 


Mg I] 


0.070 








0.08 


0.08 


27 








4571 


Mg I] 


0.068 








0.05 


0.05 


33 








4639+42 


Oil 


0.051 








0.06 


0.06 


30 








4649+51 


Oil 


0.049 








0.07 


0.07 


28 








4658 


[Fc III] 


0.047 


0.36 


0.36 


23 


0.29 


0.30 


14 








4701 


[Fe III] 


0.037 


0.10 


0.11 


43 


0.10 


0.10 


24 








4711 


[Ar IV]+HeI 


0.034 


0.44 


0.53 


21 


0.50 


0.59 


11 
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Table 3 — Continued 



Position 1 Position 2 Position 3 



A (A) 


Id. 


/(A) 


F(\) 


7(A) 


% err. 


F(\) 


/(A) 


% err. 


F(X) 


7(A) 


% err 


4740 


[Ar 1 V J 


0.028 








0.06 


0.06 


31 








4755 


[T7 TTTi 

[be 111] 


0.024 








0.07 


0.07 


29 








4861 


H/3 


0.000 


100.00 


100.00 


2 


100.00 


100.00 


1 


100.00 


100.00 


3 


4881 


n?„ TTTi 
[be 111J 


-0.004 


0.24 


0.24 


28 


0.10 


0.10 


23 








4922 


u„ T 
He 1 


-0.013 


1.10 


1.14 


13 


1.03 


1.07 


7 


1.29 


1.38 


24 


4yoy 


rn TTTi 

[U 111J 


n noi 
-U.UZ1 


1 1 P. Oft 
110. ZD 


11/1 AO 

114. 4z 


o 

z 


1 q c; on 


1 QQ QO 

loo.oy 


i 


1 0*7 1 *7 


1 O/i or: 


o 
z 


4986 


rT7„ TTTi 
[re 111J 


-0.027 


0.60 


0.59 


18 


0.30 


0.30 


i 1 

14 








5007 


m TTTi 
[O 111] 


-0.032 


345.85 


338.55 


1 


404.17 


oar c o 

396.63 


1 


375.64 


367.27 


2 


5016 


u„ T 
He 1 


-0.034 


2.48 


2.49 


9 


2.44 


2.45 


5 


2.74 


2.78 


17 


5041 


q TTT 
a 111 


-0.010 








0.11 


0.11 


22 


0.72 


0.70 


33 


5048 


xi T 
He 1 


-0.041 


0.27 


0.32 


26 


0.20 


0.25 


17 








5056 


o TTT 
a 111 


-0.043 








0.08 


0.07 


27 








5192 


r a TTTi 
[Ar 111] 


-0.072 








0.08 


0.08 


26 








5198 


rivT Ti 
[N 1] 


-0.074 


0.21 


0.20 


30 


0.12 


0.12 


21 








5270 


TT7„ TTTi 
[le 111J 


-0.089 


0.17 


0.17 


33 


0.16 


0.16 


19 








5517 


rm TTTi 
[CI 111] 


-0.140 


0.41 


0.38 


22 


0.43 


0.41 


11 


0.46 


0.43 


41 


5537 


rm TTTi 
[Ol 111) 


-0.144 


0.33 


0.30 


24 


0.31 


0.29 


13 


0.35 


0.33 


47 


5755 


rivr TTi 
[N 11] 


-0.191 


0.22 


0.20 


30 


0.15 


0.14 


19 








5876 


TT T 

He 1 


-0.216 


11.76 


10.58 


4 


12.25 


11.17 


2 


11.74 


10.82 


8 


6300 


[r~\ Ti 

[O 1] 


-0.285 








1.11 


0.98 


7 


1.99 


1.77 


20 


6312 


ro TTTi 
[S 111] 


-0.286 


1.12 


1.00 


19 


2.02 


1.78 


5 


1.12 


1.00 


26 


6364 


[O 1] 


-0.294 








0.38 


0.33 


12 


2.96 


2.62 


16 


6371 


fn TTTi 

[S ill] 


-0.295 








0.05 


0.05 


33 








6548 


nvr TTi 
[N 11] 


-0.320 


5.42 


4.61 


6 


2.94 


2.55 


4 


2.68 


2.35 


17 


6563 


Ha 


-0.322 


339.15 


288.06 


1 


330.81 


287.11 


1 


327.01 


287.21 


2 


6583 


[N II] 


-0.324 


12.06 


10.83 


4 


8.19 


7.10 


3 


8.04 


7.05 


10 


6678 


He 1 


-0.337 


3.45 


2.90 


7 


3.73 


3.23 


4 


3.81 


3.37 


14 


6716 


[SII] 


-0.342 


12.85 


10.83 


4 


10.04 


8.63 


2 


10.34 


9.00 


9 


6731 


[SII] 


-0.343 


9.46 


7.97 


5 


8.35 


7.18 


3 


7.20 


6.26 


10 


7065 


He I 


-0.383 


2.54 


2.11 


9 


2.95 


2.51 


4 


2.32 


2.02 


18 


7136 


[Ar III] 


-0.391 


8.87 


7.30 


5 


9.90 


8.34 


2 


8.86 


7.57 


9 


7281 


He I 


-0.406 


0.71 


0.59 


16 


0.75 


0.63 


9 








7319 


[OH] 


-0.410 


3.65 


2.97 


7 


3.79 


3.16 


4 


2.87 


2.43 


16 


7330 


[OH] 


-0.411 


2.98 


2.43 


8 


3.03 


2.53 


4 


2.45 


2.08 


18 


7442 


Nil 


-0.422 


0.24 


0.20 


28 


0.25 


0.21 


15 








7510 


ClI 


-0.429 


0.82 


0.67 


15 


1.06 


0.88 


7 


1.03 


0.87 


27 


7552 


[Fe II] 


-0.440 


0.81 


0.65 


15 


0.88 


0.72 


8 


0.93 


0.78 


29 


7720 


[Fe II] 


-0.450 


0.86 


0.68 


15 


1.14 


0.94 


7 








7751 


[Ar III] 


-0.452 


2.02 


1.62 


10 


2.05 


1.68 


5 


0.99 


0.83 


28 


7994 


Nil 


-0.473 


1.35 


1.07 


12 


1.38 


1.12 


6 








8360 


P22 


-0.503 








0.14 


0.12 


20 








8375 


P21 


-0.504 








0.11 


0.08 


23 








8413 


P19 


-0.507 








0.22 


0.18 


16 








8467 


P17 


-0.512 








0.36 


0.29 


12 
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Table 3 — Continued 











Position 1 






Position 2 






Position 3 




A (A) 


Id. 


/(A) 


F(\) 


/(A) % 


err. 


F{\) 


/(A) % 


err. 


F(X) 


/(A) % 


err. 


8487 


[Fe II] 


-0.514 








1.91 


1.53 


5 


2.60 


2.12 


17 


8502 


P16+[C III] 


-0.515 








0.51 


0.41 


10 








8546 


P15+[C1 III] 


-0.518 


0.43 


0.33 


21 


0.45 


0.36 


11 


1.57 


1.28 


22 


8599 


P 14 


-0.523 


0.34 


0.26 


24 


0.58 


0.46 


10 








8665 


P 13 


-0.528 


2.59 


1.99 


9 


0.80 


0.64 


8 


1.43 


1.16 


23 


8694 


N 11+ [Nell] 


-0.530 








0.13 


0.10 


21 








8751 


P12 


-0.535 








1.01 


0.80 


7 








EW(H/3), in A 






341 






383 






227 




EW aba (Kp), in A 






2.4±0.1 






2.0±0.1 






2.0±0.1 




C(H/3), 


in dex 






0.16±0.03 






0.19±0.03 






0.21±0.03 




F(H0), 


in erg s~ 1 cm -2 






1.66xl0" 13 






1.04xl0" 13 






2.46X10- 14 




HHP), 


in erg s -1 cm~ 2 






2.71X10" 13 






1.60X10" 13 






3.62X10- 14 
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Table 4. NCG 460: Line Intensities of Position 1 



A (A) 


Id. 


/(A) 


F(X) 


/(A) 


% err. 


3695 


H17 


0.264 


0.84 


1.06 


20 


3704 


H16 


0.262 


2.03 


2.42 


13 


3712 


H15 


0.260 


1.68 


2.07 


14 


3722 


H14+[S III] 


0.257 


2.35 


3.05 


12 


3726 


[O II] 


0.246 


99.14 


111.49 


2 


3729 


[0 II] 


0.255 


118.27 


132.94 


2 


3734 


H13 


0.254 


1.82 


2.52 


14 


3750 


H12 


0.250 


2.30 


2.92 


12 


3770 


Hll 


0.245 


4.74 


5.65 


9 


3798 


H10 


0.238 


5.25 


6.28 


8 


3820 


He I 


0.233 


2.66 


3.04 


11 


3836 


H9 


0.229 


6.29 


7.49 


7 


3867 


[Ne III] 


0.222 


19.41 


21.47 


4 


3889 


H 8+He I 


0.218 


17.53 


20.04 


5 


3967 


[Ne III] 


0.201 


21.89 


22.68 


4 


4026 


He I 


0.193 


2.02 


2.11 


15 


4067 


[S II] 


0.182 


1.92 


0.76 


23 


4102 


H<5 


0.176 


24.47 


26.53 


4 


4340 


H7 


0.128 


45.48 


48.77 


3 


4363 


[O III] 


0.122 


3.43 


3.62 


10 


4471 


He I 


0.094 


3.97 


4.27 


9 


4563 


Mg I] 


0.070 


0.69 


0.71 


22 


4571 


Mg I] 


0.068 


0.31 


0.32 


33 


4861 


H/9 


0.000 


100.00 


100.00 


2 


4881 


[Fc III] 


-0.004 


0.82 


0.82 


21 


4922 


He II 


-0.013 


1.21 


1.28 


17 


4959 


[O III] 


-0.021 


100.79 


99.07 


2 


5007 


[O III] 


-0.032 


303.60 


296.83 


1 


5016 


He 1 


-0.034 


2.28 


2.31 


12 


5876 


He I 


-0.216 


11.50 


10.39 


6 


6300 


[O I] 


-0.285 


2.31 


2.00 


12 


6312 


[s HI] 


-0.286 


2.20 


1.90 


13 


6548 


[N II] 


-0.320 


5.32 


4.53 


8 


6563 


Ha 


-0.322 


359.21 


305.67 


1 


6583 


[N II] 


-0.324 


15.00 


12.73 


5 


6678 


He I 


-0.337 


3.81 


3.26 


10 


6716 


[s II] 


-0.342 


27.65 


23.26 


4 


6731 


[s II] 


-0.343 


19.15 


16.10 


4 


7065 


He I 


-0.383 


2.65 


2.22 


11 


7136 


[Ar III] 


-0.391 


9.53 


7.83 


6 


7281 


He I 


-0.406 


1.37 


1.14 


16 


7319 


[OH] 


-0.410 


4.75 


3.87 


9 


7330 


[OH] 


-0.411 


3.85 


3.13 


10 


7751 


[Ar III] 


-0.452 


2.11 


1.68 


13 


EW(U/3), in A 
EW abs (H/3), in A 






285 
2.0±0.1 
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Table 4 — Continued 



A (A) 


Id. 


/(A) 


F(A) /(A) % err. 


C(H/3), 


in dex 




0.21±0.03 


^(H/3), 


in erg f 


i _1 cm 


3.97xl(r 14 




in erg s 


_1 cm -2 


6.21X1CT 14 
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Table 5. Atomic Data for CELs in IRAF v2.15.1 



Ion 


Transition Probabilities 


Collisional Strengths 


N° 


Wiese et al. (19961 


DoDita et al. ( 19761 


N+ 


Wiese et al. (19961 


Lennon & Burke ( 19941 


0° 


Wiese et al. ( 19961 


Bhatia k Kastner ( 19951 


0+ 


Wiese et al. ( 19961 


McLaughlin & Bell (19931 


0++ 


Wiese et al. ( 19961 


Lennon & Burke (19941 


Ne++ 


Kaufman & Sugar ( 19861 


Butler & ZeiDDen (19941 


S+ 


Keenan et al. ( 19931 


Ramsbottom et al. (19961 


S++ 


Kaufman & Sugar ( 19861 


Galavis et al. (19951 


C1++ 


Kaufman & Sugar ( 19861 


Butler & ZeiDDen (19891 


Ar++ 


Kaufman & Sugar ( 19861 


Galavis et al. (19951 


Ar +3 


Kaufman & Sugar ( 19861 


ZeiDDen et al. (19871 
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Table 6. Densities and Temperatures 







NGC 456 




NGC 460 




Position 1 


Position 2 


Position 3 


Position 1 


Densities (cm 3 ) 










[0 II] 


80±15 


130±30 


60±30 


170±20 


[S II] 


< 250 


250±60 


< 500 


< 150 


[CI III] 


< 3000 


< 1000 


< 1300 




Temperatures (K) 










[0 II] 


13500±250 


12400±400 


12600±800 


12600±600 


[N II] 


11300±850 


11700±1000 






[S II] 


8500±700 


9700±450 


7300±1400 


8500±450 


[0 III] 


12650±200 


12165±160 


12300±550 


12400±450 
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Table 7. Ionic Abundance Determinations from Recombination Lines 3, 



Ion 


NGC 456 




Position 1 


Position 2 


Position 3 


i z = 0.000 t 2 = 0.035 ± 0.032 


t 2 = 0.000 i 2 = 0.067 ± 0.013 


i 2 = 0.000 t 2 = 0.040 ± 0.040 


He+ 10.910±0.011 10.906+0.012 


10.923+0.008 10.918+0.008 


10.945+0.015 10.924+0.019 


C++ 


7.47+0.18 7.46+0.18 




o++ 


8.20+0.19 8.20+0.19 





NGC 460 Position 1 
t 2 = 0.000 t 2 = 0.036 ± 0.027 
He+ 10.935+0.014 10.921+0.015 
0++ 7.86+0.30 7.86+0.30 



a In units of 12+ Log N(X +i )/ N(H), gaseous content only. 
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Table 8. Ionic Abundance Determinations from Collisionally Excited Lines a 



Ion 








NGC 456 










Position 1 




Position 2 




Position 3 




t 2 = 0.000 


i 2 = 0.035 ±0.032 


t 2 = 0.000 


i 2 = 0.067 ±0.013 


i 2 = 0.000 


i 2 = 0.040 ± 0.040 


C++ 






6.75±0.15 


7.28±0.18 






N° 






4.54±0.08 


4.71±0.06 






N+ 


6.03±0.03 


6.15±0.14 


5.92±0.04 


6.07±0.04 


5.90±0.07 


6.02±0.16 


o° 






6.00±0.06 


6.15±0.05 


6.18±0.11 


6.30±0.18 


o+ 


7.31±0.03 


7.46±0.18 


7.35±0.05 


7.53±0.04 


7.38±0.10 


7.53±0.19 


o++ 


7.76±0.02 


7.87±0.13 


7.88±0.01 


8.13±0.07 


7.83±0.06 


7.96±0.17 


Ne++ 


6.96±0.03 


7.09±0.14 


7.11±0.02 


7.39±0.07 


7.10±0.07 


7.25±0.19 


S+ 


5.32±0.03 


5.44±0.14 


5.33±0.03 


5.48±0.04 


5.29±0.07 


5.40±0.15 


S++ 


6.24±0.05 


6.37±0.14 


6.30±0.03 


6.58±0.07 


6.03±0.12 


6.18±0.22 


C1++ 


4.34±0.09 


4.45±0.18 


4.40±0.05 


4.64±0.08 


4.43±0.16 


4.56±0.29 


Ar++ 


5.58±0.03 


5.68±0.11 


5.64±0.02 


5.86±0.06 


5.42±0.06 


5.54±0.16 


Ar+ 3 






3.88±0.12 


4.14±0.14 








NGC 


460 Position 1 












t 2 = 0.000 


t 2 = 0.036 ± 0.027 










N+ 


6.16±0.05 


6.28±0.12 










o+ 


7.57±0.06 


7.73±0.16 










o++ 


7.72±0.05 


7.84±0.12 










Ne++ 


7.03±0.06 


7.16±0.14 










S+ 


5.71±0.04 


5.83±0.12 










S++ 


6.30±0.08 


6.43±0.14 










Ar++ 


5.62±0.04 


5.72±0.11 











a In units of 12+ Log N(X +i )/N (H), gaseous content only. 
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Table 9. NGC 456 Position 2 Gaseous Abundance Determinations 1 



Element 


This paper 


pTp b 


PEFW C 


GISFHP d 




t 2 = 0.000 t 2 = 0.067 ±0.013 


t 2 = 0.055 


t 2 = 0.000 


t 2 = 0.000 



He e 


10.923±0.008 


10.918±0.008 


10.898 








7.57 ±0.18 


7.56 ±0.18 








C f 


6.85 ±0.10 


7.38 ±0.18 








N f 


6.67 ±0.02 


6.89 ±0.03 


6.61 


6.48 


6.43±0.02 


O c 


8.29 ±0.19 


8.29 ±0.19 








o f 


7.99 ±0.02 


8.23 ±0.05 


8.12 


8.07 


8.06±0.01 


Ne f 


7.25 ±0.02 


7.51 ±0.05 


7.22 


7.12 


7.22±0.02 


S f 


6.45 ±0.04 


6.72 ±0.13 






6.46±0.02 


Cl f 


4.53 ±0.05 


4.78 ±0.11 






4.66±0.04 


Ar f 


5.76 ±0.04 


5.98 ±0.13 






5.77±0.01 



Tn units of 12+ Log N{X)/N(R). 



^Peimbert Torres-Peimbertl (| 19 761 ). 



Pagel et al.1 (jl978l ). 



faseva et al.1 (|201lh . 
e Recombination lines. 
f Collisionally excited lines. 
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Table 10. NGC 456: Positions 1 and 3 Gaseous Abundance Determinations 3, 



Element 


Posii 


tion 1 


Posii 


tion 3 


t 2 = 0.000 t 


2 = 0.035 ±0.032 


t 2 = 0.000 t 


2 = 0.040 ± 0.040 


He b 


10.910±0.011 


10.906±0.012 


10.945±0.015 


10.924±0.019 


N c 


6.70 ±0.02 


6.79 ±0.07 


6.57 ±0.05 


6.68 ±0.09 


O c 


7.89 ±0.02 


8.01 ±0.10 


7.96 ±0.05 


8.10 ±0.13 


Ne c 


7.13 ±0.02 


7.26 ±0.09 


7.26 ±0.05 


7.41 ±0.13 


S c 


6.39 ±0.04 


6.52 ±0.17 


6.20 ±0.07 


6.35 ±0.21 


Cl c 


4.39 ±0.03 


4.50 ±0.08 


4.48 ±0.05 


4.61 ±0.08 




5.71 ±0.05 


5.81 ±0.16 


5.55 ±0.06 


5.67 ±0.19 



a In units of 12+ Log N(X)/N(R). 
b Recombination lines. 
c Collisionally excited lines. 



-36 - 



Table 11. NGC 460 Position 1 Gaseous Abundance Determinations' 



Element 


This paper 


PEFW b 




t 2 = 0.000 t 2 = 0.036 ± 0.027 


t 2 = 0.000 



He c 


10.935±0.014 


10.921±0.015 




N d 


6.60 ±0.04 


6.70 ±0.07 




O c 


8.09 ±0.30 


8.10 ±0.30 




o d 


7.96 ±0.04 


8.09 ±0.09 


8.07 


Ne d 


7.30 ±0.03 


7.45 ±0.08 


7.36 


S d 


6.50 ±0.06 


6.63 ±0.16 




Ar d 


5.85 ±0.08 


5.95 ±0.14 





a In units of 12+ Log N(X)/N(H). 



Page! et all (119781 '). 



c Recombination lines. 
d Collisionally excited lines. 
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Table 12. NGC 456, NGC 6822 V, NGC 346, 30 Doradus, the Orion Nebula, and 



Protosolar Total Abundances Relative to O a 


Element 


NGC 456 b 


NGC 6822 V c 


NGC 346 d 


30 Doradus c 


Orion f 


Suns 


12 + log He/H 


10.918±0.008 


10.909±0.011 


10.900±0.003 


10.928±0.003 


10.988±0.003 


10.98±0.01 


12 + log O/H 


8.33±0.05 


8.45±0.06 


8.21±0.06 


8.62±0.05 


8.77±0.03 


8.73±0.05 


log C/O 


-0.83±0.37 


-0.34±0.13 


-0.93±0.08 


-0.48±0.05 


-0.25±0.04 


-0.26±0.15 


log N/O 


-1.44±0.09 


-1.40±0.17 


-1.40±0.15 


-1.27±0.08 


-1.04±0.10 


-0.86±0.17 


log Ne/O 


-0.82±0.11 


-0.82±0.09 


-0.89±0.06 


-0.79±0.06 


-0.72±0.08 


-0.76±0.14 


log S/O 


-1.51±0.19 


-1.65±0.09 


-1.65±0.12 


-1.63±0.10 


-1.55±0.05 


-1.57±0.13 


log Cl/O 


-3.55±0.17 


-3.74±0.10 




-3.70±0.12 


-3.44±0.05 




log Ar/O 


-2.35±0.18 


-2.39±0.08 


-2.39±0.10 


-2.36±0.10 


-2.15±0.06 


-2.29±0.13 


log Fe/O 




-1.44±0.10 


-1.47±0.10 




-1.27±0.20 


-1.19±0.11 



a The O and C gaseous abundances have been corrected for the fractions of these elements trapped in dust grains, 
see text. 

°Values in this column are relative to O with t 2 = 0.067 ± 0.013. 

c Nebular abundances, values for t 2 = 0.076 ± 0.018, obt ained for NGC 68 22,from jPeimbert et al.ll2005l) . with the 
exception of the Fe/O value that comes from stellar data llVenn et al.ll200lTl . 

feufour et al.Nl982h ; |Peimbert et al.1 ll2000hfeelano et al.1 d2002h ; Peimbert, Peimbert. fc Luridianal lb002T) , value s 
for t 2 = 0.022±0.008. The Fe/O value comes from stellar data llVennlll999l;lRolleston et al.ll2003l;lHunter et al.ll2005h . 

^Peimbert] (120031 1. values for t 2 = 0.033±0.005. 

lEsteban et al. N2004I) . values for t 2 = 0.022±0.002. The O and C abundances have been increased by 0.08 dex and 
0.10 dex respectively to take into account the fractions of these elements trapped in dust grains. The CI abundance 
has been decreased by 0.13 dex due to an error of +1.00 dex in the determination of the Cl + /H" 1 " ratio. 

g Taken from the protosolar abundances of lAsplund et aL I l l2009h . 
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Table 13. NGC 456 and NGC 460 Total Abundances Relative to O a 



Element 


NGC 456 Position 2 


NGC 460 Position 1 


t 2 = 0.000 


t 2 = 0.067 ±0.013 


t 2 = 0.000 


t 2 = 0.036 ±0.027 


12 + log He/H 


10.923±0.008 


10.918±0.008 


10.935±0.014 


10.921±0.015 


12 + log O/H 


8.09±0.02 


8.33±0.05 


8.06±0.04 


8.19±0.09 


log C/O b 


-0.82±0.37 


-0.83±0.37 






log C/O c 


-1.24±0.12 


-0.95±0.23 






log N/O 


-1.42±0.03 


-1.44±0.09 


-1.46±0.07 


-1.49±0.17 


log Nc/O 


-0.83±0.03 


-0.82±0.11 


-0.76±0.07 


-0.74±0.17 


log S/O 


-1.64±0.06 


-1.51±0.19 


-1.56±0.10 


-1.56±0.26 


log Cl/O 


-3.56±0.06 


-3.55±0.17 






log Ar/O 


-2.32±0.06 


-2.35±0.18 


-2.20±0.13 


-2.24±0.23 



a The O and C gaseous abundances have been corrected for the fractions of these elements 
trapped in dust grains, see text. 

b Abundancc from RLs relative to O from RLs. 

c Abundancc from CELs relative to O from CELs. 
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(b) 



Fig. 1.— (la) VLT image of NGC 460 Position 1 and Positions 1, 2, and 3 of NGC 456. 
This image of 306" x 87" has Position 2 located at a = l ft 13 m 50.8 s and 6 = -73°18'03.2" 
(J2000.0), and was taken with a filter that suppresses light bluer than 4350 A. (lb) Same 
image but with a different saturation level so that Position 2 in NGC 456 is enhanced. 
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Fig. 2.— Spectra of Positions 1, 2, and 3 in NGC 456, and Position 1 of NGC 460. These 
spectra include the blue (3600-5500 A) and red (5500-7400A) high-resolution spectra, as well 
as the low- resolution spectrum (7400-8800 A). 
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Fig. 3. — Zoom-in of the region between A4550 and A4680 in NCG 456 Position 2 to show 
the quality of the data. 



